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Effects of the helicity on the dynamics of turbulent flows are investigated. The aim is to disentangle 
the role of helicity in fixing the direction, the intensity, and the fluctuations of the energy transfer 
across the inertial range of scales. We introduce an external parameter a that controls the mismatch 
between the number of positive and negative helically polarized Fourier modes. We present direct 
numerical simulations of Navier-Stokes equations from the fully symmetrical case, a = 0, to the 
fully asymmetrical case, a = 1, when only helical modes of one sign survive. We found a singular 
dependency of the direction of the energy cascade on a , measuring a positive forward flux as soon 
as only a few modes with different helical polarities are present. Small-scale fluctuations are also 
strongly sensitive to the degree of mode-reduction, leading to a vanishing intermittency already 
for values of a ~ 0.1. If the analysis is restricted to sets of modes with the same helicity sign, 
intermittency is vanishing for the modes belonging to the minority set, and it is close to that 
measured on the original Navier-Stokes equations for the other set. 


The direction of the energy transfer in a turbulent 
flow is believed to be determined by the combined ef¬ 
fects of all inviscid invariants which depend on the em¬ 
bedding dimensionality and/or on the coupling with ex¬ 
ternal fields, as in conducting or buoyant systems HHS|. 
For fully homogeneous and isotropic turbulence (HIT) 
in two dimensions, the presence of two positive-definite 
invariants, energy and enstrophy, leads to a split regime 
with energy flowing towards large scales (inverse cascade) 
and enstrophy to small scales (forward cascade) [2H3- 
Three-dimensional (3D) Navier-Stokes equations (NSEs) 
possess two inviscid invariants; energy and helicity, the 
scalar product of velocity and vorticity [101112] . Different 
from the energy, helicity is observed to be preserved by 
some dissipative events, such as antiparallel vortex re¬ 
connection ii, and it is not positive definite. As a 
result, it is not possible to predict the direction of the en¬ 
ergy and helicity transfers from fundamental arguments. 
Numerical simulations, phenomenological arguments, dy¬ 
namical models, closures, and comparison with the invis¬ 
cid Gibbs-like equilibrium distribution suggest that both 
energy and helicity develop a forward cascade in HIT 
fl2lfT5H2Tl . On the other hand, it is well known that the 
external mechanisms such as rotation [22/ (23], confine¬ 
ment (241 [25] , shear [2Bj or coupling with the magnetic 
field [27] might revert the direction of the energy cas¬ 
cade. Strikingly enough, such a reversal of the flux has 
been predicted and observed also in 3D HIT with explicit 
breaking of parity invariance, i.e., by restricting the dy¬ 
namics to a subset of Fourier modes such that the helicity 
becomes sign definite [25H50] , suggesting that inverse en¬ 
ergy transfer events are much broader than previously 
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thought and they are potentially present in all flows in 
nature. Another not understood crucial aspect of fully 
developed turbulence is intermittency , the tendency of 
the flow to develop more and more non-Gaussian velocity 
fluctuation at smaller and smaller scales. It is not known 
which key degrees of freedom are responsible for such a 
phenomenon and whether it is crucial to keep all interac¬ 
tions among different scales to preserve it. Consequently, 
we cannot model it and we cannot predict its degree of 
universality, i.e., independence of the large-scale configu¬ 
ration. There exist models based on strong mode reduc¬ 
tion (shell models [3T]) which preserve many of the in¬ 
termittent properties of the original Navier-Stokes equa¬ 
tions. On the other hand, an increase or decrease of 
intermittency is observed in numerical simulations of the 
original 3D equations with a modulation of local or non¬ 
local Fourier interactions |32] . In this Rapid Communi¬ 
cation we systematically investigate the effects of helical 
mode reduction in 3D NSEs. The aim is threefold: We 
want to explain the role played by helicity in fixing (i) the 
direction, (ii) the intensity, and (iii) the intermittency of 
the energy transfer mechanism. 

The key tool used is based on a suitable projection 
of the NSEs allowing one to disentangle, triad by triad , 
the properties of the energy transfer as a function of the 
percentage of negative helically polarized modes kept in 
the simulation. The existence of a control parameter is 
crucial to address the problem in a quantitative way, 
tailoring the degrees of freedom retained and removed, 
without any modeling. We start with the helical decom¬ 
position [16] [33] of the velocity field u(x), expanded in 
a Fourier series as uj~ = + u^h^, where h J are 

the eigenvectors of the curl, i.e., ik x = ±kh J. We 
choose hfc = 0k x k ± iOk, where 0k is a unit vector 
orthogonal to fe, satisfying the condition 0k = —0_k, 
e.g., i>k = z x k/\\z x fe||, with any arbitrary vector z. 
In terms of such an exact decomposition of each Fourier 
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mode, the total energy, E = f d 3 x \u(x )\ 2 , and the total 
helicity, H = f d 3 xu(x) ■ u>(x), are written as 

E = J2\ u i\ 2 +\ u k\ 2 ’ H = l%| 2 ), (1) 

k k 

where u> is the vorticity (see also Refs. mm for other 
previous applications of the same decomposition). The 
nonlinear term of the NSE can be then decomposed in 
terms of the helical content of the complex amplitudes, 
u s k with Sfc = ± (see Ref. UBj). We consider the dynam¬ 
ics of an incompressible flow (V • u = 0) determined by 
the decimated NSE in which a fraction a of the negative 
helical modes has been switched off [35] . We introduce 
the projector on positive/negative helical modes as 


RUN 

1-8 

9-13 

14-19 

20 

21 

22 

N 

256 

256 

512 

512 

1024 

1024 

kf 

[1.3] 

[1.3] 

[1,2] 

[42,50] 

[10,12] 

[1,2] 

a 

0 - 0.999 

0.1 -0.9 

0 - 0.9999 

1.0 

1.0 

0 

Fo 

1.0 

1.0 

1.0 

10.0 

1.0 

1.0 


TABLE I: N: Number of collocation points along each axis. 
kf. forced wavenumbers. Fo: Forcing amplitude. RUN 1-8: 
Decimation of only negative helical modes with probability 
in the range at [0: 0.9999]. RUN 9-13: Same of RUN 1-8 
but with either positive or negative helical modes (with 50% 
probability) removed. RUN 14-19: Similar to RUN 1-8 at 
higher resolutions. RUN 20-21: Forced at small scales. RUN 
22: Same as RUN 1 at higher resolution. 




h k < 8 ) h k 


( 2 ) 


where • denotes the complex conjugate. We define an 
operator D a that projects each wavenumber with a prob¬ 
ability 0 < a < 1 , 


u a {x) = D a u(x) = ^ e lkx V k Uk, (3) 

k 


where T> k = (1 — j k ) + j k F k and the random numbers 
are j k = 1 with probability a or 7 ^ = 0 with probability 
1—a. The a-decimated Navier-Stokes equations (a-NSE) 
are 


d t u a = D a {-u a -Vu a - Vp a ] + vAu a , (4) 

where v is the viscosity and p is the pressure. Notice that 
the nonlinear terms on the right-hand side of Q are fur¬ 
ther projected by D a in order to enforce the dynamics on 
the selected set of modes for all times. Despite the fact 
that the a-NSE break the Lagrangian properties of the 
nonlinear terms [35], both energy, E = J2k(\ u k I 2 + (1 — 
7 fc)|Mfc| 2 ), and helicity, H = E fe k (\ u k I 2 “ ( X ~ 7 fc)l% | 2 ), 
are still invariants in the inviscid limit of ([4]), as in any 
Galerkin truncation. We can then identify two extreme 
cases: When a = 0, we recover the original NSE, and 
when a = 1 , helicity becomes a coercive quantity with 
a definite sign. It has been recently shown jS3J[25] that, 
in the latter case, the dynamics of Q develops a dou¬ 
ble cascade characterized by an inverse energy transfer 
with a Kolmogorov spectrum E(k) ~ k~ 5 ^ 3 for wavenum¬ 
bers smaller than the forcing scale, k -C kf, and a di¬ 
rect helicity cascade with a fc ~ 7 / 3 spectrum for k kf. 
Here we address the questions: Does there exist a critical 
value a c where the direction of the mean energy transfer 
suddenly reverses as observed in two-dimensional hydro- 
magnetic systems when changing the forcing mechanisms 
E3, or does the helicity play a singular role? Are a 
few modes with opposite helical sign enough to transfer 
energy to small scales (a c —> 1), as suggested in Ref. 


)30| from considerations based on absolute equilibrium? 
What happens to small-scale intermittency in the for¬ 
ward cascade regime? Does it depend on the amount of 
negative/positive helical inodes retained? Is the residual 
small-scale vorticity mainly helical? In order to answer 
these key questions, we performed a series of numerical 
simulations at changing a with a fully dealiased, pseu- 
dospectral code at resolution up to 1024 3 on a triply pe¬ 
riodic cubic domain of size L = 2n. The flow is sustained 
by a random Gaussian forcing with 

( fi(k,t)fj(q,t')) = F(k)S{k - q)6(t - t')Qij(k), 

where Qij(k) is a projector assuring incompressibility 
and F{k ) = F 0 k~ 3 . F 0 is nonzero only for |fc| S [fc m in : 
fcmax]- See Table. [I] for details of the simulations. A dif¬ 
ferent forcing process based on a second-order Ornstein- 
Uhlenbeck process has also been implemented for some 
simulations with the same results (not shown). In all 
cases we have used a fully helical forcing with projection 
only on h k in order to ensure a maximal helicity injec¬ 
tion rate h, independent of the degree of decimation a of 
the negative helical modes. 

Energy transfer. We start from the spectral properties 
of the system following Ref. m- We define the total 
spectra restricted to the positive/negative helical modes 
as E+(k) = E|fc|= fe K\ 2 ’ E ~( k ) = E| fc |=fc( 1 -7fe)l%| 2 , 

and the corresponding quantity for the helicity, H ± (k) = 
fcE ± (fc). In the case when both energy and helicity are 
transferred forward with a rate e and h , respectively, we 
expect the usual Kolmogorov 1941 scaling (K41) for both 
energy and helicity spectra [nun], 

E{k) ~ C E £ 2/3 k~ 5/3 , H(k) ~ C H h£- 1/3 k~ 5/3 , 

which reflects in the scaling for each component as 

E ± (k) = £ 2/3 k~ 5/3 [ 1 ± Ch{£k)~\ (5) 

where C = Ch/Ce■ In Fig.[l]we show the time evolution 
of the total energy E starting from a null configuration 
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FIG. 1: Evolution of energy at varying a. All simulations 
reach a stationary state, except for the case at a = 0.999, 
where a constant increase of energy is observed, signaling the 
existence of a stable inverse energy transfer. 

at t = 0 when varying the degree of decimation from 
a = 0, for the nondecimated NS case, to cr ~ 1. We 
notice first that the time needed to develop the initial 
release of energy becomes longer with increasing a, and 
that the oscillations around the stationary regime are 
also larger when a ~ 1. The most striking phenomenon 
is that even for a very high decimation of negative he¬ 
lical modes, a ~ 1, the system is able to reach a sta¬ 
tionary state by transferring energy to the small scales. 
In other words, it is enough to have very few negative 
helical modes to develop a stable and stationary positive 
energy flux. This is quantified in Fig. [2j where we sepa¬ 
rately plot the spectra for the two helical components for 
various a. The spectrum for the positive helical modes 
[Fig. § a)] is almost unchanged and independent of a 
with a clearly developed k ~ 5 / 3 slope, whereas the spec¬ 
trum for the negative helical modes [Fig. [2^b)] tends to 
react back and become more and more energetic as a in¬ 
creases; this can be explained by looking at the behavior 
of the energy flux. In Fig. [2jc) we show that the energy 
flux is constant and independent of a for all a < 1—it 
reverts only for a ~ 1. The surprising efficiency of the 
nonlinear transfer suggests that helicity plays a singular 
role in turbulence: A tiny mixture of positive and nega¬ 
tive helical modes axccross all scales is enough to sustain 
a forward energy. This can be reconducted to the role of 
the triads with two high-wavenumber modes of opposite 
helicity nu. If this is the case, the most important triads 
must have one negative and two positive helical modes: 

S(k\p,q) = ((k-u-)(u+-u+)} + ({k-u+)(u+-u~)). ( 6 ) 

They are present with probability oc (1 — a) while triads 
with two negative helical modes, exist with probability oc 
(1—a) 2 . To keep the above triadic correlation constant at 
decreasing a , we must have: uj~ —> u jT/(1 —a) and there¬ 
fore: E~(k) = E|fc|=fc( 1 ^ lk)Wl | 2 E~(k)/(l - a). 

This prediction is shown to be well realized in the inset 




FIG. 2: (a) Log-log plot of E + (k) = X[|fc|=fc l M fc | 2 vs k at 
changing a. (b) Log-log plot of E~ (k) = J2\ k \ =k ( 1 -‘lk)\u^\ 2 
vs k at changing a. Inset: Rescaled E~ (k) with factor (1 —a), 
(c) Semilog plot of flux of energy. Inset: Flux of helicity, at 
changing a. 


of Fig. gb), where we show that rescaling E (k) by a 
factor (1 — a) leads to a good overlap, except for a ~ 1, 
where the fluctuations due to the onset of the inverse en¬ 
ergy transfer becomes very large and the above argument 
possibly breaks down. Thus, negative helical modes act 
as a bridge for the energy transfer; they receive energy 
from the large-scale positive modes and release it to the 
small-scale positive modes; the fewer there are, the more 
intense their amplitude must be to do it efficiently. When 
negative helical modes become too rare or absent, i.e., for 
a ~ 1, this bridging is no longer possible and the energy 













































4 


flows upscale [25]. Helicity plays the role of a passive cat¬ 
alyst in the energy transfer. Proving the existence of a 
unique a c for the inversion of the energy transfer could be 
extremely hard and it may not be crucial. The observed 
value is so close to unity that it might also be dependent 
on the realization of 7 *, and/or on the Reynolds numbers. 
This issue is left for a more detailed analysis in a future 
work. 



FIG. 3: (Color online) Iso-vorticity surfaces for (a) a = 0, 
(b) a = 0.5. The color palette is proportional to the intensity 
of the helicity, red for high positive values (~ 10 3 ) to blue for 
high negative values (~ — 10 3 ). 

Intermittency. The second important problem addressed 
concerns intermittency, the presence of strong non- 
Gaussian fluctuations at small scales, usually interpreted 


as a build up of instabilities in the vortex-stretching 
mechanisms. Here, we want to understand how inter¬ 
mittency changes under helical mode reduction. A vi¬ 
sual inspection of the vorticity field, in Fig. [3] shows a 
strong depletion of filament-like structures, of the stan¬ 
dard 3D NSE [Fig. [3ja)] with decimation of the nega¬ 
tive helical modes, shown, e.g., for a = 0.5 [Fig. |3][b)]. 
In order to quantitatively assess the degree of intermit¬ 
tency, we focus on the so-called structure functions (SFs), 
S'(p)(r) = ((S r u a ) p ), based on the p moments of the trans¬ 
verse velocity increments S r u a = Uy{x + r) — Uy(x) as 
a function of the separation scale (the selection of the 
x — y components is arbitrary because of isotropy). In 



FIG. 4: Top panel: Excess kurtosis of transverse SF for r = rj 
at changing a: (□): decimation of negative helical modes 
only; (0) : decimation of either positive or negative helical 
modes with 50% probability; (A): a posteriori decimation of 
negative helical modes from a velocity held of standard non- 
decimated NSEs. Bottom: Logarithmic local slopes for fourth 
order moment, a = 0, i.e., the full NSE (□). a = 0.5 : The 
whole held (O) or its h £ (A) and h ^ (0) components. The 
horizontal line represents the standard intermittent correction 
as predicted from Ref. f.’Ts . 

Fig. a we show (i) the local slopes of the relative scal¬ 
ing of fourth order SF with respect to the second or¬ 
der, C (4) ( r ) = dlog[S^ 4 \r)]/dlog[S^ 2 \r)] for the case 
a = 0.5, and (ii) the value of the excess kurtosis for 
r = 7], K(rj) = (rf)/[S^- 2 \r ))] 2 , at changing a. In 

the top panel, concerning the kurtosis, we found that 
intermittency is highly sensitive to a decimation; it is 
enough to remove a small fraction of negative helical 
modes from the dynamics to strongly deplete the non- 
Gaussian character. In the same panel we show also the 
results of another numerical experiment, where we re¬ 
peated the measurements in a set of simulations (RUN 
9-13) with random decimation; this time either a positive 
or a negative helical mode is decimated with a probabil- 
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ity a. The reduction in the intensity of intermittency is 
comparable with the previous case. To further investi¬ 
gate the role of dynamic helical mode-reduction, we per¬ 
formed a projection aposteriori, applying the operator 
D a to the velocity field obtained from a fully resolved 
non-decimated NSE (a = 0). In this case, intermittency 
remains almost unchanged, independently of a, suggest¬ 
ing that only the dynamical mode reduction is crucial to 
deplete the vortex-stretching mechanism. In the bottom 
panel we show the results concerning the local scaling 
exponent (r) for a = 0.5, and we compare it with the 
original NS case (a = 0) and with two other important 
measurements obtained by taking the velocity configura¬ 
tions dynamically generated with a = 0.5 and projecting 
them on their positive or negative helical components, 
i.e., by applying the projector or V on all modes. 
Doing that, we observe intermittency for the projection 
on the majority component h £ and a vanishingly small 
correction for the projection on the minority helical com¬ 
ponent h ^. In summary, we can conclude that intermit¬ 
tency as measured from real-space velocity configurations 
is the results of a highly non-trivial and entangled cor¬ 
relation among a subset of key modes in Fourier space. 
It is hidden in the correlation among all Fourier modes 
with given helical components. 

Conclusion. We have highlighted and quantified the sin¬ 
gular role played by helical Fourier modes in the energy 
flux reversal, showing that a forward transfer is always 
preferred as soon as a very small percentage of modes 
with opposite helicity are present. These findings sug¬ 
gest the possibility to check a posteriori on direct nu¬ 
merical simulations and experiments of strongly rotating 
flows, of flows under vertical confinement or under strong 
shear the role played by different helical triadic interac¬ 
tions in driving the energy transfer forward or backward. 
Intermittency measured on real-space fields is fragile and 
strongly dependent on the mode-reduction protocol, sug¬ 
gesting that its origins must rely on highly nontrivial cor¬ 
relations among helical and nonhelical fluctuations. In 
particular, it is apparently key to have all modes with a 
given helical component. Another key factor might be 
to keep local or nonlocal interactions in Fourier space, 
as also suggested by numerical experiments based on a 
scale-dependent mode-reduction scheme [52: . 
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